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Abstract 



O 

(R-channel) TBA is elaborated to find the effective central charge dependence on the 
boundary parameters for the massless boundary sine-Gordon model with the coupling 
constant (8n) / (3 2 = 1 + A with A a positive integer. Numerical analysis of the massless 
boundary TBA demonstrates that at an appropriate boundary parameter range (cusp 
■ point) there exists a singularity crossing phenomena and this effect should be included in 
TBA to have the right behavior of the effective central charge. 

>< 

1 Introduction 

The low dimensional quantum system such as a quantum wire with boundaries is not 
easy to study in terms of mean field approach due to large quantum fluctuations. The 
system is also strongly affected by the the existence of boundaries. For example, one 
needs a good knowledge of the the low dimensional quantum field theory to study the 
quantum Hall edge tunnelling pQ . 

In this work, the massless Tomonaga-Luttinger liquid with boundaries is studied mo- 
tivated by SNS junction super-conductor analysis[21 E]. This system is summarized in 
terms of the boundary sine-Gordon model(bSG). 

A= I ' d 2 x —(d a <p) 2 - $ I rfx2cos(%-^ 1) ))-/ig ) / cfe2cos(%-^l 2) )). (1-1) 
J 4it J y =o Jy=R 

The free energy dependence of the finite system on the boundary parameter x — b(Vo ~ 
(p^) is one of the demanding questions. To do this, we first consider the massive sine- 
Gordon model with boundaries and put the bulk mass vanish 
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The massive bSG is written as 



A 



d x 



2/x(cos(26(^) - ] 

H$ [ dx2 cos(% - y^)) - /4 2) / dx2 cos(% - ^ 2) )) , (1-2) 

Jy=0 Jy=R 



The coupling constant b 2 is restricted to be less than 1. (Note that b 2 is scaled by 8ii from 
the conventional choice (3 2 = 8irb 2 ). 

The bulk sine-Gordon model (SG) [S] belongs to the category of two dimensional in- 
tegrable quantum field theories and allows an exact treatment of the system. Integrable 
quantum systems have been studied systematically after the pioneering work of Zamolod- 
chikov 0. The integrability of the bSG was also demonstrated in jjj. 

The scale dependence of the system can be studied by the method of thermodynamic 
Bethe ansatz (TBA) [El El- Suppose a system lies along the y-axis with a finite size R and 
appropriate boundary conditions are imposed at each end as in Fig |T|). The x-direction 
is periodic and its size L is put to oo in the thermodynamic limit. 




Figure 1: space with two boundaries: 
boundary condition (a) and boundary 
condition {(3). 
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Figure 2: space with states: initial 
states \B a ) and final state \Bp). 



The partition function with these boundaries is given as 



Tre 



-LH, 



-LE aP {R) 



for large L. H a p is the Hamiltonian of the system of size R with boundary (a) and (/?). 
E a p(R) is the ground state energy in the thermodynamic limit, which depends on the size 
R. 

The same system can be viewed as the one with initial state \Bp) and final state \B a ). 
In this picture, the bulk is periodic in L as in Fig. (J2J). Then the same partition function 
is evaluated using the periodic Hamiltonian H of the system. 



-RLfcf){R) =<Ba \e- RH \B g >= 



< B a \A >< A\B b > e~ REA 
< A\A > 



(1-3) 
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where f a p{B) is the free energy density per length and {A} is the complete set of the 
periodic Hamiltonian eigenstates. The finite size effect of the SG (with boundary) was 
analyzed in |UEI| for diagonal case using thermodynamic Bethe ansatz (R-channel TBA). 

In section 2, we present the massive (R-channel) TBA for the bulk sine-Gordon model 
with boundary sine-Gordon interaction. In this analysis, we restrict the coupling constant 
A = 1/b 2 — 1 = rib + 1 to a positive integer (n& > 0) so that the bulk scattering matrix 
is diagonal but non-diagonal boundary scattering is allowed. The topological charge 
violation at the boundary is incorporated following the suggestion given in [2|. This 
TBA has the bulk and boundary scale dependence as well as the boundary parameter 
dependence. 

In section 3, (R-channel) TBA for the massless boundary sine-Gordon model is ob- 
tained as the massless limit of the massive TBA. It is demonstrated through numerical 
work that beyond a certain parameter range (% > b 2 ir), a branch singularity crossing 
occurs. Thus the original ground state TBA which holds within a small parameter range 
is to be modified. We investigate on the singularity structure of TBA and find that only 
one singularity of each species is moving on the complex plane. This has the effect on 
the evaluation of the convolution in TBA. From this branch singularity crossing informa- 
tion, we propose a modified massless TBA according to the suggestion in HH H2 EH] 
similarly in the bulk case [Hj. Section 4 is the conclusion. 



2 Massive TBA for boundary sine-Gordon 

The bulk sine-Gordon periodic potential allows a soliton with topological charge +1 
and an antisoliton with the charge —1. Both of them have the same mass M as the result 
of the charge conjugation symmetry, — > —(f). The mass is given in terms of \x [15J, 
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(2-1) 



In addition, there are topologically neutral particles, breathers (interpreted as the soliton- 
antisoliton bound states). Their masses are given as 

M a = Mm a = 2M sin a = 1,2, ■ ■ ■ ,n b . (2-2) 

n& is the number of breather species, rib = positive integer less than A. 

The free energy density in Eq. (jl-3|) is expanded in terms of the bulk Hamiltonian 
eigenstates, i.e., solitons, antisolitons and breathers. These states are uniquely identified 
in terms of mass and rapidity due to the Fermi-statistics. 

The presence of the boundary forces two restrictions on the states. First, the pair 
creation at the boundary forces the rapidity paired (6, —6) and therefore, one can count 
9 positive and make energy eigenvalue doubled: Ea{0) = 2MacosIi9 where Ma is the 
single particle mass. 

Second, at the boundary an in-coming soliton is allowed to be scattered away as an 
antisoliton and vice versa since the soliton number is not conserved in general. To take 
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Figure 3: pair creation at the boundary 



care of this, soliton and antisoliton are regarded as a constituent of a doublet of identical 
particles .2j. When the partition function in Eq. (jl-3|) is written in terms of spectral 
density p, Z a/3 = J[dp]exp(—RL f a/3 (L, R)), the spectral density should include not only 
one-particle density of topological particle (denoted as p ) (soliton is indistinguishable 
from antisoliton) and its hole density (pho), one-particle density of a breather (p a ), its hole 
density ( ph a with 1 < a < n ) but also should include two-particle density of topological 
particles (pd) {i.e. soliton and antisoliton pair). The densities are summarized in the 
Tabled 



species 


density 


soliton or/and antisoliton 


Po, Pd, PhO 


breather 


Pa, Pha, a = l,---,n b 



Table 1: particle species and the corresponding densities. 



Then the free energy density can be written as 
Rf a p{R)= Hdel J2 PA(2M A cosh8-\n\^(6))+ Pd (2M d cosh6-\n\ d a p)-S B 

J ° {A=0,l,-.n b 

(2-3) 

where M = M and M d = 2M. S B is the entropy density, 

S B = / d0[ (p + Pd + Pha) hi(p + Pd + Pha) ~ Po In Po - Pd In p d - pho In pm 
Jo 

+ {(Pa + Pha)^(pa + Pha) - Pa^Pa- Pha^Pha}] ■ (2-4) 

0=1, •••,71;, 

Arg(0) = (B a \9A) (6A\Bp) is the boundary state contribution, which is given in terms of 
the boundary scattering amplitude R{u) with u = —id: 



Kxp = K-{u)K b p{u) foro = l,2,...,n 

= K^K+ + + K*=K+- + K^Kp + + K^Kjf = Tr(K a K p ) 
\% = (K++ K-~ - K+- K-+)(K++Kp- - K+~Kp+) = Det(iQ^) (2-5) 

where K(u) = R(n/2 — u). 

The boundary scattering amplitude (modulo CDD-type factors), can be found in 
which satisfies the boundary version of the Yang Baxter equation, unitarity condition, 
and analyticity- crossing symmetry. 

- ( R ++ ^,u) R + ~(v,$,u) ) 



P + (v,$,u) Q (u) 
Qo(u) P (r],$,u) 



Ro(u 



a(rj, u) cr(it?, u) 



cos(?7) cosh($) 

Pq(j], u) = cos(Aii) cos(?7) cosh(i?) =p sin(A-u) sin(^) sinh($) 
Qo(u) = — sin (Aw) cos(Ait) . 

Here Rq is the boundary condition independent part, 

r(4ZA-^)r(4A(Z-l) + l 2Xu " 



(2-6) 



Ro(u) = n 



1 _ ' ^ U 



-u 



J((4/ - 3)A - ^)r((4Z - 1)A + 1 - ^) 
and cr(x,u) is the boundary condition dependence part, 

T(i + a + (2/ - 1)A - £)r(| - £ + (2Z - 1)A - £ 



COS X 



a\x,u) 



cos 



(x + Am) 2 



r(| - a + (2/ - 2)A - f )r(| + a + 2ZA - f ) 



"/(« 



-u 



The scattering parameters, 77 and $ are related with the action parameters, and (fio 



«yWsin(& 2 7r) cos(6v?o)/v / ^ 



cos(& 2 r/) cosh(& 2, #) 
sin(6 2 ^) sinh(fe 2 ^) = /i B y / sin(6 2 7r) sin(6<yj )/ v //I 

The boundary scattering amplitude of breathers is given as 

(77, 0, it) = i?f } («) (*7, ^ W «) , 1 < A; < n b . 

i?Q fe is the boundary independent part and the boundary dependent one: 

k-21-1 



(2-7) 



(2-8) 



2A 



(A + i) " (A + f) 

where the notation (x) stands for 



1=0 



(JL 4. I 4. 



2A 



sin (| + 



sin (f - f ) 



The hole and the particle densities are not independent each other. This relation is 
obtained from the bulk scattering amplitude. The bulk-scattering amplitude of solitons 
and antisolitons [Sj are given as 



Stt(u) 

siz(u) 



S (u) = s(u) 



S+-(u) 



sin(Aw) 



sin(A(7r — u)) 

sin(A7r) 
sin(A(7r — u)) 



s(u) 
s{u) 



(2-9) 
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where s(u) is given as 

T(2(/ - 1)A - ^f)T(2l\ + 1 - ^; 



s(u) = n 

1=1 



r((2/-i)A-^)r((2/-i)A + i-^ 



Due to the restriction of A, the bulk scattering amplitudes are diagonal, S'+l(w) = 0. 
This restriction makes our analysis not too much complicated j2D|- The diagonal scattering 
amplitude for soliton and antisoliton turns out to be equal up to a phase difference: 

sn(u) = (-i) x - 1 suu). 

The scattering amplitude of the breathers B a and B b with b < a < rib takes the form 
S ab (u) = {a + b- l}{a + b - 3} ... {a - b + 3}{a - b + 1} , (2-10) 
where the notation {y} is defined as 

(y±l\ (y-l\ 

and satisfies the relations {y}{— y} = 1 and {y + 2X} = {—y} . The scattering amplitude 
of the soliton (antisoliton) and breather S^(u) = S%+(u) = S^Z(u) is given as 

q (a) ( \ _ f {a-l + A}{a-3 + A}---{l + A} if a is even 

{) { -{a-l + A}{a-3 + A}---{2 + A}(4±!)(4^) if a is odd ' 1 } 

Demanding the wave function periodic in L we have the constraints between hole 
densities with particle densities. For soliton states (n , n d , Uho) we have 

N t 

exp(zLM sinh 0°) J[ S 00 (9° - fyS^ + flJ^ooM)^? - 9 d )S 0d (9° + Vf) 
II (SUA ~ 0j)SM + 0])j I = ±e 2m ^^ +n ^ +n ^) . (2-12) 
For breathers (n a , n a h, a = 1, 2, • • • , m): 

N ( 

exp(rLM a sinh 9<*) J[ S aa (0- ~ 0])S aa (0- + 9])S aa (29?)S a0 (e? - 9°)S a0 (9? + 0°) 



Sad(9i — 0j)S a d(9i + 0j) Y[ [ S a b(9i ~ 9j)S a b(0i + 9j 



6=1 



= ±e 27ri(n a (e?)+n fea (0?)) _ (2-13) 

Differentiating with respect to the rapidity, we have the relations of hole and particle 
spectral densities: 

roc 

M AC oshe+ Yl / d9 B PB(0 B )(<f>AB(e A -e B ) + <pAB(e A + e B )) 

B=0,l,-,n& 

+ / de d Pd (6 d ) (0 Ad (e A - e d ) + + o d )) 

Jo 

= 2n(p A (9) + p hA (6) + 6 A0 p d {9)) for A = 0, 1, ■ ■ • , n b , d, (2-14) 
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where p A (9) = <p AB (9) = -j ^alg) and = 2<f> AO (6). 

Introducing pseudo energies, e 

e —=P±. foro = l,...,n 6j e- eo = ^, e" 6 " = ^ , 

P/ia P/iO P/lO 

and minimizing f a p{R) we have the massive TBA: 

e A {0) = d A (o)-]n\£g-— £ / de'<t) AB {e-e')L B {e'), 

271 B=0,l,-,n b J -°° 

e d (9) = 2(e + lnA° /3 )-lnA^, (2-15) 

where d A (0) = 2 M A R cosh 6 , L = ln(l + e" e ° + e" £d ) , and L a = ln(l + e" ea ) ( a = 
1, •••,n&). 

Shifting the pseudo energy eA — > 6 a — hi Aa, we have = 2e and the boundary TBA 
is put equivalent ly as 

e A (9) = d A (9) - — E / d9> MO - 9') L B {9') (2-16) 

ZH g J— oo 

where A = 0, 1, • • • , n& and 

£o s lna + A^e^ + A^e- 2 -) 

L a = ln(l + A^e" €a ) for a = 1, • • • , n b . 

In terms of this TBA, the free energy has the form, 

7T 

E(R) = Rf(R) = -R£bulk + ^boundary ~ TTT^Ceff (2"17) 

where £buik = — i^ 2 ta n ^x is the bulk energy density, ^boundary is the boundary energy 
obtained by Al. Zamolodchikov ^H]; whose details can be found in ^7j 

. „, M ( .77. , 1 . 7T , 1 . , 7T . l\ 

boundary ( V , *) = " ^/JA) { + COSh(-) - - COs(-) + - Sm(-) - - j , 

and c e fj is the effective central charge 

6RM r°° - 

c eS {RM) = — — d9 m A cosh 9 L A (9) . 

n J ~°° A=0,l,-,n b 

One may also put this TBA in a reduced form as in a bulk TBA ^H] for rij, > 1. To do 
this, we include the doublet formally by extending the whole index into A' = 0, 1, • • • , n&, + 

e A ,(0) = d A ,(0) - — £ / d9'<t> A , B ,{9-9')L BI {9') (2-18) 

B =0,1, ■••,«.(,, + 

where + denotes the doublet, e + = e d /2, m + = m , and L + = . Then using an identity 
of the fourier transformed form of the kernel, (p A > B i{k) = J^° 00 d9 <p A > B '(9) e lke , 

(W-^^W) = 5 ^'- 2c o^/h) lA ' B ' (2 " 19) 




where h = 2 A = 2n b + 2 and Xa> b> is the incidence matrix of D nb+2 , |U HHJ d] (Ia' b' — 1 
if yl'-B' is directly connected, otherwise), 



1 2 s — 1 

o--o- -o- -o- 



we can put the TBA in a reduced form when n& > 1, 

/oo , 
^'K(# - e'){L B , + e fl > - d B ,){0') , (2-20) 

J B'=U,l,-,n i) ,+ ~°° 

where K(6) is the reduced kernel, 

K ^ = 2tt cosh(A I) 

This reduced TBA allows one to construct Y-system equations with Y{9) = e e ^ \ 



17T ITT 

Y Q (e + T )Y (e- T ) = [KM + YnM] 



Yn b {0 + -)Y nb (6 - -) = [\ d (9) + \ (9)Y (9) + Y (9) 2 ][K b ^ 1 (9)+Y nb . 1 
Y a (9 + -)Y a (9 - -) = H[\ b (9) + Y b (9)] x ^ for a = 1, 2, ■ ■ ■ , n b - 1 

Y + {9) = Y (9). (2-21) 

3 Massless TBA for boundary sine-Gordon model 

The massive TBA shows the scale dependence as well as the parametric dependence. 
Since we are particularly interested in the parametric dependence of the boundary TBA 
on 0o, we will look at the massless limit of the TBA in this paper. 

The massless TBA corresponding to bSG action Eq. (jl-lj) is obtained by taking the 
limit — > of the massive TBA in Eq. (|2-lfi[ 12-20)1 . Even though the soliton mass M 
vanishes, one may introduce a finite renormalized mass scale Mr as Mr = (M/2)e 6 ° if a 
large parameter 9 is defined as 



6 2 \ 2A/(A+1) 

-2b 2 ) 

7(* ! ) \V^n^. 



(. 



In this limit, the rapidity is rescaled into renormalized one, 9r as 9 = Or + 0q. The 
boundary scattering parameter, d is also rescaled as $r maintaining a relation •§ — X9 = 
$r — X9r. Then $r is written in terms of the action parameters, 



m R f x e b2 *« = 2/WCo sin(&%) . (3-1) 
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On the other hand, rj is not rescaled but is identified as 

h 2 n = i b< ^° fOT - ~ 71 'I 2 ' C\ ?\ 

' \ 7T - b(j)o for 7T/2 < 600 < 7T . 1 1 

This identification is justified from the numerical analysis later on. In terms of this 
rescaled parameter (we omit hereafter the subscript R standing for the renormalized one) 
the solitonic boundary scattering amplitudes are given as 

K h ^ u ^{ _ ie ~W) ^a,) J ^ X/4 *(«) 0-3) 

where u = n/2-u and kiu)" 1 = UqLo 2 cos (§ + g + , and the breather bound- 

ary reflection amplitudes are given as 

K W (tf , ,,„)_» TT sin(|(^-f -f-vr^l^)) 
The boundary \ a p is given as 

yd ^ _ | e (i? Q -A6»— i7rA/2) _|_ e -(i? a -Ae-t7rA/2) j, 

jg^-A^A/2) + j (^kju) k/3(u)y (3-5) 

From this massless scattering data the massless TBA is given as the form 

1 roc 

e A (9) = d a (o) - — / de'4>AB(e-e')L B {e') (3-6) 

Z7T g J —00 

where D A > = 2m A /re e with r = MR and the soliton mass M in M^/ is replaced by 
the renormalized mass Mr. The reduced form of TBA can be similarly obtained. TBA 
Eq. (|3-6|) . however, turns out to be better suited for the numerical analysis of various 
range of parameters. The free energy is given as 

71 

Rf{R) — ^boundary — 24^> Ceff (3 _ 7) 

where the bulk energy density vanishes, and the boundary energy and effective central 
charge are given as 

boundary 2 COs(7r/2A) 



c eff = — / dO Yl ™ A e 9 L A (6) 

IT J ~°° A=0,l,-,n b 



Let us investigate the parametric dependence of the ground state energy Eq. (13-71) 
using TBA Eq. (13- For the details of analysis, we will restrict our selves to A = 1 and 
A = 2 cases with the symmetric boundaries ($ Q = dp = d). 
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When A = l(n{, = 0), the boundary contribution is given explicitly as, 



A af3 



cost] e 



( 2 



4cosh(^)cosh(^ 



V/3 



tanh( 



0, 



9, , 



tanh(- 



TBA is trivial since the kernel 



0: €o(0) = 2m§re e . The effective central charge 



is obtained numerically and is plotted c e ff v.s. x i n Fig. HI where x — Kfa 



(2) 



may put 



(One 



using the freedom of field translation). We note that c e ff is 7r-periodic 



in x- This is because the boundary fugacity A° # in Eq. (|3-5|) is 27r-periodic in r\ = rji — r] 2 
and X a/ 3 is r\ independent. Generally, c c g- will be 2nb 2 - periodic in x- 

However, the periodicity of the energy in x is n °t acceptable as pointed out in [2]. 
When x > 7r /2, the boundary term in the Lagrangian effectively changes the relative sign; 
one can equivalently put //^ — > —fig and — > //^ while X ^ ^ — X- This relative 
sign change of the boundary term should be reflected in the c e g value. The same problem 
was also observed in boundary Lee- Yang model 





Figure 4: c e g vs. x when A = 1 and d = 10 
before modifing the massless TBA. 



Figure 5: c c g vs. x when A = 1. Modified 
TBA shows the correct behavior (in stars) 
for = 10 (sold) and $ = (dashed). 



To cure this disease it has been proposed in OHHj that a certain excited state contribu- 
tion should be properly taken care of. According to this proposal, the c e ff is recalculated 
and is presented in Fig. 03 when A = 1. In the analysis two things are considered for 
X > v/2: First, one needs to find the zeroes of e L ° in terms of the complex rapidity 9 
following JT] with the parameter identification in Eq. (|3-2|) . 

1 + \% e~^ §) + \% e~ 2 ^ §) = . (3-8) 

Eq. ()3-8|) is solved numerically. It turns out that there are infinite number of solutions for 
9 = ±iir/2+9p with 9 P real and therefore, the pseudo energy is given as €q(9) = ±i2more 6p . 

Second, one needs to take account this branch singularity effect into the free energy. To 
understand which solution has the effect on the free energy for x > 7r/2, we just complexify 
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the value of % (equivalently O ) around n/2 (cusp point) as x ^ tt/2 — (^)e l7rfc with k 
varying from — > 1 and use the TBA Eq. ()3-6|) allowing 9 complex but 9' real to trace the 
movement of singularities. With the notation t = e L ° = l + \° a0 e~ e ° + X^ e~ 2e ° , the result 



is given in Fig. [H] for $ =. We choose the case of $ = since this does not correspond 
to the conformal boundary condition, neither Dirichlet not Neumann condition. (Note 
that the singularity position is contour plotted using 1/(1 + |t |) rather than l/|to| to 
renormalize the height). 
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Figure 6: Singularity crossing for A = 1 case: Contour plot of 1/(1 + |to|) is presented in 
the complex 9 plane. The vertical axis stands for — | < Im(#) < | and the horizontal 
one for < e Rc ^ < 0.5. From the top left to the bottom right, each figure is reproduced 
for X = tt/2 - ^e ink with k = 0, 1/5, 2/5, 3/5, 4/5, 1 for = 0. 

When x — real, singularities lie on the imaginary value lm(9) = ±iir/2. Only one of 
the singularities (the left-most one with the smallest value of 9 P ) crosses the real rapidity 
line during the complexified changes of x horn the negative imaginary value — in/ 2 to the 
positive imaginary plane and finally sits at the positive imaginary value m /2. (When k 
varies from to —1, the left-most singularity at m /2 is moving down and crosses the real 
axis). The other signularities remains in the imaginary axis with Im(#) = ±7r/2. This is 
seen in Fig. |3 Similar behavior is demonstrated in scaling Lee- Yang model ^TJ El E! 
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and used in boundary SLYM (Crossing singularities are paired for massive theories). 




15.5 




15.5 




15.5 




10.5 15.5 




10.5 15.5 




15.5 



Figure 7: Other singularities for A = 1 case: Contour plot of 1/(1 + |to|) is presented in 
the complex 9 plane. The vertical axis stands for — ~ < Im(#) < | and the horizontal for 
0.5 < e^ ^ < 16. From the top left to the bottom right, each figure is reproduced for 
x = tt/2 - ^e i7rk with k = 0, 1/5, 2/5, 3/5, 4/5, 1 for ■& = 0. 

Due to this singularity crossing, the rapidity integration in the free energy has to 
detour the branch singularity. As the singularity crosses the real axis and moves up 
to Im(#) = 7r/2 and one has to evaluate the contour integration around the singularity 
6 = in/ 2 + 6 P . The result is given as 



— c eff = -im e 6 - — 4 d6m e e L (9) 
24r 2tt J-oo 



(3-9) 



where represents the principal value of the integration. The the first term in the RHS is 
the singularity crossing contribution and the second one is the real axis contribution. (Here 
we adopt the convention of the singularity crossing from the negative imaginary value 
to the positive imaginary one. One may equivalently consider the opposite singularity 
crossing and detour the integration down to the negative imaginary complex plane and 
get the same result. Note that if 6 is a solution, so is 6*). 
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We generalize this idea into the case A > 1. To understand the singularity crossing, 
we need to know the singularity structure in the comlex 9 plane. Note that the TBA 
Eq. ()3-6|) is singularity- free on the rapidity domain with —ir/h < lm(8) < ir/h when \ is 
real. One can analytically continue the TBA into other domain of 9 using the relation of 
Y-system Eq. (122T1) . 

When A = 2 the boundary \ a p is given as 



2 


cos?7 e ( 2 - a > -f e 1 2 Z0 J 




[4cosh(^" 2e + cosh(^- 2e l l) 




^cosh^" 29 + cosh(\ 2e 7)] 



d 4 cosh(tf a - 29) cosh(^ - 20) 

aP ~ (4cosh(^ + f ) cosh(^ - f )) 2 (4cosh(^ + f ) cosh(^ - f )) 2 ' 
Xlp = tanh(#/2 - # a /4) tanh(#/2 - ^/4) , (3-10) 

where 77 parameter is identified according to Eq. (j3-2j) and the kernel is given as 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 



Figure 8: Singularity structure of TBA in the complex 9 plane when A = 2 and x = 
Ti / 3 — 7r/30: Contour plot of 1/(1 + \t \) (up) and 1/(1 + \zi\) (down). The vertical axis 
stands for — | < lm(9) < j and the horizontal one for < e Re ^ < 6. The t zero with 
the arrow attached is moving. 

c c ff is 7r-periodic in \ an d nas the cusp at 7r/3. (See Fig. HU below). Fig. El is the 
numerical study of the singularity in the domain with lm(0) < tt/2 at real x — vr/3 — 7r/30 
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Figure 9: Contour plot of 1/(1 + \t Q \) (up) and 1/(1 + |zi|) (down) for k = 0.25 (left) and 
k = 0.5 (right). The vertical axis stands for — f < lm(8) < f and the horizontal one for 
< e Rc ^ < 6. The moving zeros of t and Z\ are pointed with arrows. 

below the cusp point. Here t = e L °. In the figure Z\ = e ei and 1/(1 + \ t \) and 1/(1 + \zi\) 
are plotted after normalizing l/|to| and l/|zi|- 

It is to noted that all the zeros of to lie at Im(6 l )=±7r/2 instead of Im(#)=±7r/4. In 
addition, there are no zeros of z\ in this domain. (Note that pseudo energy e A is symmetric 
under 9^9* when \ is real). 

Next, the singularity positions are traced with the complexified %■ X is P u t as X ~ y 
f — (if)) e * 7rfc anc ^ ^ * s var i e d — > 1. The numerical result is given in Fig. El As increases, 
one of the t zeros (with the smallest value of real rapidity) goes up. When k = 0.25, 
one zero of z\ with the same real rapidity appears at Im(#) = —n/2 and moves up as k 
increases. As k reaches 0.5, the zero of Z\ comes to lm(9) = — 7r/4. 

In fact the positions of zeros of each species are related with the positions of zeros 
of other species according to the Y-system Eq. (|2-21|) . As the singularity of one species 
moves around, the corresponding singularities of other species follow the trace. 

The position of to zero is drawn in Figs. ()10| fTTjl . The data is taken up to k = 0.5 
without any difficulty However, beyond k = 0.5 numerical work becomes harder due to 
the numerical instability. We use the method using the "damped" iterative method used 
in JI] to easy this instability up to near k = 0.56. Beyond that point numerical instability 
does not allow any further computation. Neverthelss, the data clearly indicates that the 
imaginary position changes linearly in k. 

Now, it should be noted that as the to zero crosses the real rapidity line, we cannot use 
the TBA Eq. (J3-6J) any longer. This is because the moving singularity pushes up the real 
integration line of the convolution and forces to modify the TBA. Denoting the position 
of the to zero as 9, we have 

~ 1 r°° ~ 

e A (9) = 2m A re e + In S A0 (9 -9)-— ]T f d9'<p AB (9 - 9')L B (9') . (3-11) 

Z7T B=0,l J -°° 

Here the branch-cut contribution is used as follows. 

__L [ de'<j> A0 (6 - 9>)L(9') - \nS M {9 - 9) - -L / d9'<p M (9 - 9')L(9') . 
2,71 J 2ix J 
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Figure 10: t crossing singularity position: 
e Re ^ v.s. k. Data (denoted as o) for less 
than k = 0.6 are taken from the original 
TBA. Data (denoted as *) for close to k — 
1 are taken from the modified TBA. 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Figure 11: to crossing singularity position: 
Im (9) v.s. k. Data (denoted as o) for less 
than k = 0.6 are taken from the original 
TBA. Data (denoted as *) for close to k — 
1 are taken from the modified TBA. 



Using this modified TBA Eq. ()3-ll|) we obtain Figs. (fT2l IT^j) for the singularity struc- 
ture. The figures show that to zero and z\ zero have the same real rapidity [e^ cd « 0.1) 
but with the imaginary rapidity tt/A apart. In addition, the position of to zero in Figs. (|10l 
[TUJ) (denoted as *) ends up at the imaginary rapidity 7r/4. 
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Figure 12: Position of t zero (up) and Z\ 
(down) for k = 0.9. Vertical axis stands 
for 0.58- f < Im(6») < 0.58 and horizontal 
one for < e Rc( ^ < 1. The moving zeros 
of to and z\ are pointed with arrows. 



Figure 13: Position of t zero (up) and z\ 
(down) for k = 0.975. Vertical axis stands 
for 0.73 — j < Im(#) < 0.73 and horizontal 
one for < e Rc< ^ < 1. The moving zeros 
of to and Z\ are pointed with arrows. 



As k approaches 1, \ becomes real and the to zero reaches Im(#) = 7r/4. At the 
same time the zero of Z\ reaches the real rapidity axis, Im(#) = 0. Therefore, at real 
X = vr/3 + 7r/30 > 7r/3, the z\ zero in the convolution integration modifies TBA again. 
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For real 9 the pseudo energy becomes 



2m A re e + \nS A0 (9 - 



lnS A1 ( 



roo 



2n 



i de'^ABie-e^LBie 1 ) 

> n i J— CO 



5=0,1 



2m a re H — In 



Sac 



1 roo _ 

2 / d9'<i) AB {e-e')L B {e') 



2~S AO (0-e*) 2n B tZ,J-°° 
In the last identity the bootstrap relation of the scattering matrix is used: 
S A0 (6 - iir/A)S M {9 + in/ A) = S A1 (6) , A = 1,2. 

Explicitly, 

'9 + i7r/2 



(3-12) 



Sqq{9) = i tan 



1 ' 1 ; 1 + V2cosh(0 + i7r/2) 



Note that the convolution integration in Eq. ()3-12|) refers to the principal value due 
to the singularity of Li(6') at 9 = 9 P . The numerical integration of the convolution can 
be done by shifting the contour integration. Another way of avoiding this numerical 
instability is to use TBA with the reduced kernel. Using the kernel relation (j2-19j) . we 
have 

e A >(6) = D A ,(9)+P A ,(9)+ T a>b> -P d9' K(9-9')(L B ',+e B ,-D B ,-P B )(9') (3-13) 

B'=0,1,+ J -°° V 7 

where P A >(9) = \ In and D A , = 2m A ,re 6 '. Note that (h x + e 1 )(9') vanish at 

9 = 9„. 




Figure 14: c e g vs. \ when A = 2 and $ = 0: Solid line is obtained from the original TBA 
and dotted one from the modified TBA. 



Now, the crossing singularities contribute to the free energy: 

— 4- d9 Y, ra A e 9 L A {9) 

J— CO a n i 



7r m n 
Cog = % e e * 



24r 



V2 



(3-14) 



A=0,1 
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where 9 V is the real part of 9. This singularity position should satisfy the constraint 
equation, 

t {6 p ±i~)=z 1 (6 p ) = Q (3-15) 

Thus obtained c e g is compared with the original result in Fig. El The smooth behavior 
of c e ff in Fig. El indicates that the modified TBA is in the right track. 

Finally, we remark that the modified TBA (J3-12) I3-I3j) is also written in another form 
if the kernel relation (J2-19)) is used for the S-matrix: 

e A ,(9) = D A ,(9)+ £ T A , B , [Z(9,9) + f°° d9' K{9 -9'){Lb, + e B , - D B ,){6') \ (3-16) 

B'=0,1,+ V J -°° ) 

where 

i 5(0,0) = tan _1 (tanh(3(0 - 6)/ 2)) - tan -1 (tanh(3(0 -9 + itt/2)/2)) . 
This TBA satisfies the same Y-system in Eq. (12-21)1 with the fugacities replaced with 



the massless one (h=4 when A = 2). Redefining Y's, yo(9) = Y Q (9)/y X d (9) , yi(9) = 
Yi(9)/X 1 (9) and uo = in/ 4, we may put into this form. 

y (9 + u)y (e-u) = l + yi (6) 

yi (9 + u) yi (9-u) = l + -±= yo (9)+(y (8)y (3-17) 



where fugacity relations for A = 2 are also used, 



\ d a p(9 + tu)\ d ap (9-u) = (a^(0)) 
\lp(9 + u) \i p (9 - u) = Xt p (9). 



2 



Note that since A°(^) / J\ d {9) is — ^-periodic in 9, so are yo(9) and yi{9) as in the 



2 

bulk case. 



4 Conclusion 

We presented the (R-channel) TBA for the (massless) boundary sine-Gordon theory 
with coupling parameter A = positive integer, which incorporates the violation of the 
topological charge. With this TBA, the boundary parameter effect on the effective central 
charge is investigated. 

This TBA, however, doest not reflect the boundary parameter condition faithfully 
beyond the cusp point x — b 2 n. We have demonstrated that this unpleasant feature is due 
to the branch singularity crossing in the TBA equation. Numerical analysis show that near 
the cusp point of c e fr some of the singularities of TBA crosses the real rapidity. It turns 
out that only one singularity for each species changes its position in the complex plane. 
Especially, its imaginary position moves "linearly" with the phase of the complexified 
boundary parameter \ near the cusp. We expect this singularity crossing feature will 
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be very universal. When A = 2 (h — 4 and periodicity | in unit of 7r), the imaginary 
parts of the original singularity positions are ±|, | . After singularity crossing some of 
the positions change into 0, ±j. This leads into consideration for other A's (h = 2A). 
From the original singularity positions ±f , ±|, - • • , ±£, ^±1 with periodicity ^±2, some 
of the singularities will presumably rearrange into 0, ±i ±|, • • • , This expectation 

needs to be checked further as well as the behavior for various parameter range. The scale 
dependence of the massive TBA is also not well understood. This issues will be carried 
on in a separate paper. 

This analysis is also expected to hold for an arbitrary coupling constant case. However, 
the TBA of massless/massive boundary sine-Gordon theory with arbitrary coupling is not 
feasible at this moment since the bulk Hamiltonian eigenstates will result in the infinitely 
coupled TBA equations [20]. Instead of this approach, DDV type equation is expected 
to be more suitable, which does not impose the string hypothesis for the structure of the 
roots of Bethe ansatz [21] . 

It is known [22] that the sign change effect of the boundary term of boundary Liou- 
ville theory and boundary sinh-Gordon model [2H] can be explained using the analytically 
continued boundary parameter only to a certain range. Nevertheless, the boundary La- 
grangian phase difference is also known [22] to induce the branch singularity crossing 
similarly as in this boundary sine-Gordon model. Furthermore, the integrable boundary 
ADE-affme Toda theories [21] have discrete boundary conditions, (+), (— ) and Neumann 
condition. Among the three, (— ) boundary condition is not yet fully understood. It re- 
mains to be seen if there is any relation between the relative sign change of the boundary 
term \i B — > — ji B and the exicted state contribution to c e fj. 

Note added: While our manuscript being revised, a preprint by Caux et. al. ( |cond-m at /0306328 ) 
appeared where similar singularity structure information at the cusp point is used for the 
analysis of the finite size effect. 

Acknowledgement 

We thank C. Ahn, P. Dorey, K. Moon and R. Tateo for invaluable discussions and 
KIAS for hospitality. CR thanks Al. Zamolodchikov for explaining his work on the singu- 
larity crossing in boundary sinh-Gordon model, H. Saleur for informing their singularity 
conjecture of TBA and Durham University where the revision of the manuscript was 
made. This work is supported in part by the Basic Research Program of the Korea Sci- 
ence and Engineering Foundation Grant number R01-1999-00018-0(2002) and by Korea 
Research Foundation 2002-070-C00025. 

References 

[1] P. Fendley, A. Ludwig, and H. Saleur, Phys. Rev. Lett. 74 (1995) 3005; Phys. Rev. 
B 52 (1995) 8934. 

[2] J.-S. Caux, H. Saleur and F. Siano, Phys. Rev. Lett. 88 (2002) 106402. 



18 



I. Affleck, J.-S. Caux and A. M. Zagoskin, Phys. Rev. B 62 (2000) 1433. 

P. Fendley, H. Saleur, and N. P. Warner, Nucl. Phys. B430 (1994) 577. 

A. B. Zamolodchikov and Al. B. Zamolodchikov, Ann. Phy. 120 (1979) 253; 
L. Fadeev, E. Sklyanin and L.Takhtajan, Th. Math. Phys. 40 (1979) 688; V. E. Ko- 
repin, Th. Math. Phys. 41 (1979) 953. 

A. B. Zamolodchikov, Int. J. Mod. Phys. A4 (1989) 4235. 

S. Ghoshal and A. B. Zamolodchikov, Int. J. Mod. Phys. A9 (1994) 3841; S. Ghoshal, 
Int. J. Mod. Phys. A9 (1994) 4801. 

Al. B. Zamolodchikov, Nucl. Phys. B342 (1990) 695. 

C. N. Yang and C. P. Yang, J. Math. Phys. 10 (1969) 1115. 

A. LeClair, G. Mussardo, H. Saleur, and S. Skorik, Nucl. Phys. B453 (1995) 581. 
P. Dorey and R. Tateo, Nucl. Phys. B482 (1996) 639. 
P. Dorey and R. Tateo, Nucl. Phys. B515 (1998) 575. 

P. Dorey, A. Pocklington, R. Tateo and G. Watts, Nucl. Phys. B525 (1998) 641. 

V. V. Bazhanov, S. L. Lukyanov, A. B. Zamolodchikov, Nucl. Phys. B489 (1997) 
487. 

Al. B. Zamolodchikov, Int. J. Mod. Phys. A10 (1995) 1125. 
Al. Zamolodchikov in 5th Bologna Workshop (2001). 
Z. Bajnok, L. Palla, G. Takacs, Nucl. Phys. B614 (2001) 405. 
Al. B. Zamolodchikov, Phys. Lett. B253 (1991) 391. 

T. R. Klassen and E. Melzer, Nucl. Phy. B338 (1990) 485; B350 (1991) 635. 

M. Takahashi and M. Suzuki, Prog. Th. Phys. 48 (1972) 2187. 

C. Destri and H. J. de Vega, Phys. Rev. Lett. 69 (1992) 2313; Nucl. Phys. B438 
(1995) 413. 

Al. Zamolodchikov, in private communication. 

A. B. Zamolodchikov and Al. B. Zamolodchikov, Nucl.Phys. B477 (1996) 577-605; 
V. A. Fateev, A. Zamolodchikov, and Al. Zamolodchikov, "Boundary Liouville Field 
Theory I. Boundary State and Boundary Two-point Function", hep-th/0001012 

[24] E. Corrigan, P. Dorey, R. Rietdijk, and R. Sasaki, Phys. Lett. B333 (1994) 83-91; 
V. Fateev, "Normalization Factors, Reflection Amplitudes and Integrable Systems" , 
|hep-th/0 1030141 Mod. Phys. Lett, bf A16 (2001) 1201; C. Ahn, C. Kim, C. Rim, 
Nucl. Phys. B628 (2002) 186. 



19 



